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New antenna path-length error and pointing-error structure optimization codes 
were recently added to the MSC/NASTRAN structural analysis computer program. 
Path-length and pointing errors are important measures of structure-related an- 
tenna performance. The path-length and pointing errors are treated as scalar dis- 
placements for statics loading cases. These scalar displacements can be subject to 
constraint during the optimization process. Path-length and pointing-error calcula- 
tions supplement the other optimization and sensitivity capabilities of NASTRAN. 

The analysis and design functions were implemented as “DMAP ALTERs ” to the 
Design Optimization (SOL 200) Solution Sequence of MSC-NASTRAN, Version 
67.5. 

f. Introduction 

For many years, the principal antenna structural design tool of the Ground Antennas and Facilities 
Section has been the JPL-IDEAS Design Optimization Program [1]. It has proved very useful both 
as a design tool and as a test-bed for numerous analytical and optimization techniques. It has also 
taken advantage of numerous unique features of the FORTRAN-V compilers of the UNISYS machines. 
While this was expedient, it has proved increasingly restrictive. Costs have increased, maintainability 
has suffered, and the code is not “portable.” 

JPL-IDEAS depends on a limited repertoire of finite elements, both for computational efficiency and 
to effect the optimization algorithms. The NASTRAN structural analysis computer program has tra- 
ditionally been used to check IDEAS results and for analyses that required a more extensive element 
library. Having to maintain parallel finite-element models was at best cumbersome and time consuming. 
With the advent of design optimization capabilities in MSC/NASTRAN, and the rumored imminent de- 
parture of the UNISYS platforms, it was decided to implement the essential features of IDEAS directlv 
into MSC/NASTRAN. 

II. Analysis Capabilities 

A. Surface Accuracy Calculations 

Antenna primary surface accuracy can be represented as the weighted root-mean-square (rms) error of 
a paraboloid best fitting the deformed reflector surface in a half-path-length sense [2,3]. A subset of the 
finite-element model grid points is chosen to represent the primary-reflector panel support points. Each 
of these grid points is assigned a weighting function proportional to the reflector panel area it supports. 
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The following is a summary of the path-length error calculations for one loading case. Matrix and 
vector quantities are noted with their dimensions in square brackets, where N is the size of the weighted 
grid subset. Quantities designated as “raw” refer to errors with respect to the undeflected, nominal 
paraboloid; those designated as “fit” are errors with respect to the best-fitting paraboloid. The other 
symbols are defined as 


/ = nominal focal length of main reflector 
X> Y ? Z = spatial coordinates of the main reflector support points [N x 1] 
U, V, W = translation displacement components [ N x 1] 

wgt = diagonal matrix of surface panel weighting function [N x N] 
'Ewgt = sum of the surface panel weights 

p = raw path-length errors, with respect to nominal paraboloid [ N x 
gX, gY, gZ — direction cosines of surface normal vectors [AT x 1] 

A = basis functions of normal equations, wgtAS = wgtp [N x 6] 

S — solution of normal equations [6 x 1] 

<fr — fit path- length errors, with respect to best-fit paraboloid [N x 1] 
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The six rows of the solution vector S are the X , Y, and Z translational offsets of the vertex of the 
best-fitting paraboloid; the change in focal length of the best-fitting paraboloid; and the rotational offsets 
of the best-fitting paraboloid about the x- and y-axes. 


The coefficient matrix ([A T wgtA]), the right-hand side (wgtp), and solution (S) are printed and 
other results of the path-length analysis are tabulated for each loading case. 


B. Pointing Error Calculations 

The pointing error is calculated as the sum of the error contributions from each of the optical compo- 
nents of the antenna: the primary reflector, the secondary reflector (S/R), and the microwave-feed phase 
center [4], The errors due to the displacements of the primary reflector are calculated from the fitting 
parameters of the best-fit paraboloid (BFP). The motions of the secondary reflector are taken to be those 
of a grid point whose motions approximate the motions of the secondary. Similarly, the motions of the 
feed phase center are taken from the motions of a designated “feed” grid point. The spatial coordinates 
of these grids are not considered during pointing-error calculations. The symbols are defined as 


K = 
K 2 = 
M = 
A,C = 
F = 

dX, dY , dZ = 
dF = 
rX,rY — 


beam deviation angle of primary reflector 
beam deviation angle of secondary reflector 
magnification factor 

Cassegrain geometry parameters, where hyperboloid eccentricity = C/A 
focal length of nominal paraboloid 
translational offsets of BFP 

change in focal length of main reflector paraboloid (/bfp — f nominal) 
rotational offsets of BFP or secondary reflector 


Angular errors due to primary (BFP) motion in a symmetric Cassegrain system are calculated as 
follows: 


Motion component 

X-axis error 

Y-axis error 

Lateral translation 

-K/F X dY 

K/F x dX 

Rotation 

(1 + K)x rX 

(1 + if) X rY 


Errors due to secondary reflector motion in a symmetric Cassegrain system are calculated as follows: 


Motion component 

X-axis error 

Y-axis error 

Lateral translation 
Rotation 

(K/F - K2/F/M ) X dY 
-(2 x (C - A) x K/F ) X rX 

— (K/F - K2/F/M ) x dX 
— (2 x (C - A) x K/F) x rY 


Errors due to feed phase center motion in a symmetric Cassegrain system are calculated as follows: 
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Motion component 

X-axis error 

Y-axis error 

Lateral translation 

{K/M - 1/F) x dY 

-{K/M - 1/F) x dX 


No pointing error accrues from axial motions or focal-length change of the components in a symmetric 
system. 


The pointing error about the elevation axis can be predicted relative to the elevation- axis position 
transducer. The rotation of the grid point, whose rotation represents that of the transducer, is subtracted 
from the total pointing error about the x-axis and is reported as the “corrected” pointing error. 

The default antenna-pointing analysis calculations assume a symmetric Cassegrain configuration: a 
paraboloidal main reflector and a hyperboloidal secondary reflector with coincident foci, and a point 
source microwave feed located at the other focus of the hyperboloid. Simple descriptions of this geometry 
are sufficient input for this traditional DSN case. For non-Cassegrain or asymmetric (e.g., clear-aperture) 
antenna configurations, the pointing-error calculations shown above can be overridden for each of the 
motions of each of the optical components. Components whose override factors have values of 0.0 (the 
default) will be calculated as outlined above. If the factor is other than 0.0, the pointing error will be the 
factor times the specified displacement. 

Pointing errors about each axis are tabulated for each loading case with the results of the path-length 
analysis. f ’ 


III. New Input Data 

Although many NASTRAN input data (over those required for a simple statics run) are required to 
characterize the antenna for the new program, the amount of additional work required of the user has 
been kept to a minimum by automating the data preparation. Four simple data generation programs are 
available, each of which is written in standard FORTRAN. 

The identification of the main reflector surface panel-support grids and of grids that represent the 
motions of the secondary reflector, feed phase center, and elevation transducer is accomplished through 
the use of NASTRAN DMIG, USET, and PARAM bulk data cards [5]. 

Two other programs can be used to generate the “DESVAR” and “D VP REIT” NASTRAN cards, 
required for sensitivity and optimization, from a list of element property cards. One of the programs 
recognizes PROD, PBAR, PBEAM, and PSHELL cards and generates cards whose design variable-to- 
property relation is cross-sectional area for the bar elements and thickness for the shells. The second 
program recognizes PBAR and PBEAM cards and generates design variables and relations suitable for 
rectangular tubes. The outside dimension and wall thickness are the design variables. These programs 
can be adapted easily to process other properties or output other design variables. 

These additional data cards need only be generated once for a particular antenna configuration. The 
small amount of additional input (6 through 16 PARAMeters describing antenna geometry and pointing 
coefficients and four data cards per design constraint) can be generated easily by hand. 


Surface accuracy and pointing error responses are each calculated as scalar values for each design 
oading case. A NASTRAN scalar point is designated for each design constraint and used as a place- 
holder for these responses. Displacement design constraints are defined on the displacements of these 
scalar points. 
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IV. Program Flow 

The new antenna analysis and sensitivity code have been implemented into MSC/NASTRAN, Version 
67.5. The new code has been installed and tested on the JPL Cray Y-MP supercomputer and a Silicon 
Graphics (SGI) Crimson model workstation. 

The general flow of the NASTRAN design optimization solution sequence (for statics) is 

(1) Read and collate input data. 

(2) Update geometric and structural matrices and tables. 

(3) Solve for static responses. 

(4) Determine sensitivity of objective and constraint functions to design variables. 

(5) Perform optimization. 

The objective is usually to minimize structure weight, with constraints on antenna performance and 
element stresses. Steps (2) through (5) are iterated until the design problem converges. 

The new antenna analysis and design code are inserted into the NASTRAN design optimization so- 
lution sequence in two places: at the end of statics analysis, step (3) above, and after the sensitivity 
analysis calculations, step (4) above. 

At the end of step (3), the new code performs some input data checking, initializes and saves invariant 
antenna analysis data, and echoes various antenna description and processing parameters to the print 
file. At this point, all statics response quantities are available and static antenna analyses are performed. 
The resulting scalar quantities are inserted into the appropriate displacement vectors for each loading 
case. 

Design sensitivity analysis [6], step (4) above, computes the rate of change in structural response to 
changes in design variables. Small perturbations to the design variables are approximated in the structural 
matrices, and response derivatives are explicitly solved for. The result, for the statics case, is a set of 
displacement derivative vectors for each design variable for each loading case. Path-length and pointing- 
error analyses are performed on each of these displacement vectors, as appropriate. The original statics 
errors are subtracted to give the change in antenna performance for variations in the design variable. 
These resulting scalar quantities are inserted into the displacement vectors for each loading case and 
passed to the optimizer. 

The added antenna analysis and design codes do not modify or interfere with the operation of the 
original optimization code and are compatible with all the other sensitivity and design options in the 
design optimization solution sequence. The advent of geometry optimization in the next release of 
MSC/NASTRAN is especially intriguing. The method employed is extensible to constraints on other 
performance and operational measures, e.g., secondary reflector offsets or geometric theory of diffraction 
(GTD) performance analyses. 


V. Effectiveness 

JPL-ANTOPT is intended eventually to supplant the JPL-IDEAS program. Test design optimization 
runs have been made using a small verification model (Fig. 1) and the IDEAS design model for the new 
DSS-24 34-m antenna (Fig. 2). These models are summarized in Table 1. (The numbers in square brackets 
are the number of design variable linking groups in each element type. Path-length error constraints (P/L) 
are in mm rms. Pointing errors (PNT) are in arc-sec.) 
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Table 1. Comparison model summary. 


Component 

Test model 

34-m model 

Grids 

28 

1473 

Elements 

Rod 

88 [12] 

3876 [243] 

Tri 


106 [10] 

Quad 


562 29 

Mass 

18 

1127 

Design constraints 

^-gravity P/L 

0.254 

0.381 

T-gravity P/L 

0.254 

0.381 

y -gravity PNT 

50.0 

900.0 

Rim +Z force PNT 

50.0 


S/R +y force PNT 

100.0 



Design constraint values were chosen to exercise the programs and have no other significance. For the 
baseline configurations, each constraint is violated by a factor of from 1.5 to 4. Stress and other side 
constraints were disabled. The default optimizer input parameters were used for both programs except 
that termination of design iterations due to meeting convergence criteria was prevented. 

The design variables for each model are the cross-sectional areas of the rods and thicknesses of the 
plates. Minimum and maximum design variable limits were imposed: for the test model rod areas, 
64.5 mm 2 and 6451.6 mm 2 , and for the 34-m model, 322.6 mm 2 and 32258.0 mm 2 for the rod areas and 
2.54 mm and 254.0 mm for the plate thicknesses. Some of the structure model was excluded from the 
optimization problem, including lumped masses that were used to account for the weight of the main 
reflector panels, secondary reflector, and counterweight; the remainder is deemed the variable structure. 

Three comparative design optimization runs were performed using each model. The baseline configu- 
ration for the test model was uniform rod areas of 645.2 mm 2 ; that of the 34-m model was 3225.8-mm 2 
rod areas and 25.4-mm plate thicknesses. In addition, starting points of minimum and maximum design 
variable sizes were used. The geometry and element connectivity were not changed. 

Figures 3 and 4 are plots of the variable and feasible structure weights as functions of design iteration 
history for the six optimization runs for each program. Feasible structure weight is the variable structure 
weight scaled by the largest ratio of performance measure to design constraint. For example, if a design 
constraint were violated by 10 percent at a particular iteration, the feasible weight would be 110 percent 
of the variable structure weight. The ratio of feasible to variable weights, then, is proportional to the 
worst constraint violation. The individual curves in the plots are designated I-WGT and I-Feas for 
variable structure weight and feasible weight as reported by the IDEAS program, and A-WGT and 
A-Feas for the weights reported by the ANTOPT program. (The differences between the initial feasible 
weights of the IDEAS and ANTOPT 34-m model runs is attributable to the differences in the weight and 
stiffness characteristics of the plate elements in the two programs.) It can be seen that the new program 
consistently achieves a feasible design in a few iterations. 

JPL- ANTOPT comprises over 1800 lines of MSC/NASTRAN DMAP source code and has been ac- 
cepted by NASA’s Computer Software Management and Information Center (COSMIC) as program 


VI. Summary and Conclusion 

The JPL- ANTOPT program performs microwave antenna performance analysis and structure opti- 
mization in an efficient manner that is easy for the structure designer to use. The new code is portable 


289 




14 


12 


- 10 


X 

CD 

LU 

5 

LU 

X 

ID 

H 

O 

D 

X 

L— 

03 



6 



0 1 1 I I I | L 

INITIAL 1 2 3 4 5 6 

ITERATION HISTORY 


Fig. 3. Test model weight history: (a) initial = baseline; (b) initial = minimum; and 

(c) initial = maximum. 
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across computer platforms and is compatible with other NASTRAN optimization capabilities. The 
method employed to introduce antenna performance measures into the NASTRAN optimization sequence 
is easily adaptable for use with other performance measures. 
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Errata 

Due to a TDA Progress Report software problem, the plots in Figs. 1, 3, 5, 6, 7, and 11 of S. R. Tyler’s 
“A Trajectory Preprocessor for Antenna Pointing” that appeared in The Telecommunications and Data 
Acquisition Progress Report 42-118, vol April-June 1994, pp. 139-159, August 15, 1994, were distorted. 
Also, some numbers in the x-axis of Fig. 11 were transposed. The correct figures are provided below. 



Fig. 1. Matching of antenna and target positions for 
constant accelerations (elevation only). 



Fig. 3. Matching of antenna and target positions for 
elevation only (chasing target). 



Y A9 v o 
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Fig. 5. Area of valid solutions (a m = a max/2 and 

a o = a f= 0 ). 
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